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Based on the formalism of thermo-field dynamics a new approach for 
studying collective excitations in hot finite Fermi systems is presented. Two 
approximations going beyond the thermal RPA namely renormalized thermal 
RPA and thermal second RPA are formulated. 

1. Introduction 

A standard technique for treating quantum many-body systems at finite temperature 
T is the thermal Green function (Matsubara) method. We here outline an alternative 
approach, based on the concepts of thermo-field dynamics (TFD) [I], @], to treat 
collective excitations of hot finite Fermi systems. In this context TDF has two 
appealing features: 

- temperature effects arise explicitly as T-dependent vertices, providing a good 
starting point for various approximations 

- a generalization to the time-dependent situation is easy, since temperature and 
time are independent variables in TFD (real-time formalism). 

As will be detailed below, both features allow for straightforward extensions of 
well-established zero-temperature approximations. Some aspects of the TFD appli- 
cation to hot nuclei have been discussed previously by Hatsuda and Tanabe ||. 

2. TFD: Basic Elements 

The extension of quantum field theory to finite temperature requires a doubling of 
the field degrees of freedom. In TFD this doubling is achieved by introducing an 
additional tilde space [IJ 0. A tilde conjugate operator A is assigned an operator A 
(acting in ordinary space) through the tilde conjugation rules 

(AB)~ = AB; {aA + bB)~ = a* A + b*B , 

where A and B represent ordinary operators and a , b are c-numbers. The aster- 
isk denotes the complex conjugate. The tilde operation commutes with hermitian 
conjugation and any tilde and non-tilde operators are assumed to commute or an- 
ticommute with each other. The thermal Hamiltonian 7i = H — H is the time- 
translation operator and the excitation spectrum of the system is obtained by the 
diagonalization of 7i. 



The temperature-dependent vacuum |\I/ (T)) is the eigenvector of H with eigen- 
value 0: 

h\mt)) = q 

If one determines the thermal vacuum state as 



where E n , \n) and \n) are eigenvalues, eigenvectors and their tilde counterparts of 
the Hamiltonian H, respectively, the expectation value (^/q(T)\0\^/q(T)) will exactly 
correspond to the grand canonical ensemble average <C O ^> of a given observable O. 



3. Temperature-dependent Fock space in HF 

Restricting ourselves to two-body interactions, the iV-body Hamiltonian is written as 

H = Y1 ti2 a i a 2 + t V 123 4afa^a 4 a 3 , (1) 

12 4 1234 

where a + , a are fermion creation and annihilation operators. The one-body part 

contains the kinetic energy matrix T\i as well as the chemical potential A. 

Using the Wick theorem one can expand the Hamiltonian (1) in normal order 
(: ... :) with respect to the temperature-dependent ground state \^ (T)) and it has 
the following form: 

H = h + /in + h 22 , (2) 
where ^ 

ho = Y, ^12P21 + Vi234P3lP42 , 

12 ^ 1234 

hu =: X]*12 a l"°2 + Y Vl234fi420iCl3 '■ , 
12 1234 

/l 2 2 ='■ 7 X/ ^1234 a ^ a 2~ a 4 a 3 : • 
4 1234 

We denote the density matrix as 

Pij = (MT)\4ai\MT)) ■ 

To diagonalize the /in-part of H which is quadratic in a + ,a one has to perform the 
unitary transformation ||: 

a t = Y D lk a k i a l=Y °lk a k , 
k k 



where DD + = D + D = I and a + ,a are new fermionic ( "quasiparticle" ) operators. 

Requiring tin to be diagonal in a + ,a, a system of equations for the matrix 
elements D as well as the single-quasiparticle energy e\ is obtained: 

E I *12 + 2 ^1324P43 £>25 = £5-^15 (3) 
2 V 34 / 

with the auxiliary condition 

T,Pn = N - (4) 

i 

which implies conservation of the average number of particles. In the density matrix 
P12 = (MTM^IMT)) = J2 D *2kDi k n k 

k 

the quasiparticle occupation numbers n& = (^>o(T)\otloti\tyo{T))5ki can be deter- 
mined once an explicit expression of the thermal vacuum is specified. 

In practice it is impossible to find the exact thermal vacuum for the full Hamil- 
tonian (2). In setting up approximation schemes the usual starting point is thermal 
mean-field theory (the thermal Hartree-Fock approximation (THF)). In this case 
the thermal vacuum is an eigenvector of the uncorrelated thermal Hamiltonian : 

(/in - h u )\0(T)) = ^Eiiafai - 5+a,)|0(T)> = (5) 

i 

The solutions of eq. (5) define the vacuum |0(T)) for so-called thermal quasiparticles 

pj: 

Pi = XiOi - Vicxf 

Pi •^i^-i Ui^i 

A|0(T))=A|0(T)) = 

where the coefficients x iy i/i denote the thermal Fermi occupation probabilities of the 
states a+|0): 

Xi = yfl - fi ,Vi = yJ7i 



Jl l + expfo/T) 

The chemical potential A is adjusted to fulfill eq. (4) ensuring that the average 
number of particles in the system remains fixed. 

In terms of the quasiparticle operators a + , a the THF- vacuum state 

\0(T)) = H(x i + y i a + a + )\0)\0) 

i 

is similar to the BCS-vacuum for systems with pairing. 



Based on the mean field we are now able to construct a temperature-dependent 
Fock space. The Fock space is generated by acting with thermal quasiparticle cre- 
ation operators /3 + ,/3 + on the thermal vacuum 

/3 + |0(T)),/3 + |0(T)),/5 + /3 + |0(T)),... 

After the transformation to thermal quasiparticles the thermal Hamiltonian Ti takes 
the form: 

H = Tin + 7-^22 + 7^40 + 7^04 + 7^3i + (6) 
where TC mn consists of terms of the type (f3 + ) m (f3) n . 

1 

^22 = j ^1234(^1X2X3X4 - ymvsvdiPtM P^fo - MM PM + 

4 1234 

+ ^1234(2:12/22/3X4 - yi^XzVilPtM (3 2 p 4 

1234 

^0 4 = 7E ^1234(XlX 2 ?/32/4 ~ ym^l) PlfoPiPd 
4 1234 

^13 = \ £ ^1234(XlX 2 X 3 ?/4 - 2/12/22/3X4) (P£ ^^2^1 ~ M PlfoPi) 
^ 1234 

H.31 = Tl 13 , 7^40 = 7^04 

where 

£^1234 = E ^5678 ^51^62^73^84 • 
5678 

Given the Fock space the solution of any statistical problem is thus reduced to a 
diagonalization of the thermal Hamiltonian. 



4. TRPA 

Based on the above formulation, a straightforward way to solve the RPA at finite 
temperature (TRPA) is to apply the equation of motion method |J. As for T = 0, 
one has a Raleigh- Ritz variational principle with the thermal Hamiltonian (6): 



(tfoCOl \SQu, [H,Qt]] \MT)) = Vv(MT)\ \SQu,Qt] \MT)) 



(7) 



The potentially exact statistical variational problem (7) cannot be solved in 
practice and the class of variational functions has to be restricted. The simplest 
approximation amounts to the following trial wave function (phonon wave function): 



Qt = Y, r^PtM - 012M = E <2^ 2 - 0^12 



(8) 



12 



12 



If, in addition, the exact thermal vacuum in (7) is replaced by thermal HF-vacuum, 
the exact matrix equation (7) reduces to the usual TRPA equations for the ampli- 
tudes 024 an d the excitation energy w v \ 

£24^24 + ^22/4$^ C/l234 (^12/3013 + ^lV^l) 
^ 13 



"^42/2 ^1234(^32/1013 + ^l2/3^3l) = 
^ 13 



21 



£42042 + ^42/2 J! ^1234(^12/3013 + ^32/lV', 
^ 13 



3i ; 



-^^22/4 X ^1234(^32/1013 + S l2/3?A 
Z 13 



3l) 



—a; 



1^42 



At this stage only the terms Tin, H22, ^40 an d 7^04 play a r °l e an d this part of TC 
becomes diagonal in the phonon operators. 



H = "£u u QtQ„+n 



13 



H 



31 



In the TRPA the Q + ,Q operators are boson operators (the same holds for the 
bifermionic operators Af 2 and A12) The index v runs over TRPA solutions with 
positive and negative energies to v . The negative-energy solutions appear naturally 
and ensure that the set of thermal phonon states is complete within TRPA approx- 
imation: 

00 

1 = Y.Y.(Q») k \MT))(MT)\(Q») k , 

v k=0 

where \^o(T)) is the vacuum for thermal phonons, i.e. Q u \^o(T)) = 0. In terms of 
thermal quasiparticle operators, the TRPA vacuum is given by 



1 TT 1 

= ex P 7 



1234 



|0(T)> . 



5. Beyond TRPA 

The derivation of the TRPA, presented above, suggests that improvements known 
at zero temperature can be taken over straightforwardly. Below we demonstrate this 
for the so-called renormalized RPA [6-9] as well as the second RPA [TTH 



Within the renormalized RPA, the bifermionic operators Af 2 and A12 are no 
longer bosonic, as in RPA, but obey modified commutation relations: 



Am, A 



12, ^34 



^13^24 (1 - qi - 92) , 
where qi are c-numbers. At finite temperature these are defined as follows: 



The new thermal vacuum state \^/q(T)) is different from the Hartree-Fock or TRPA 
vacuum state and will be determined below. 
The substitution || [F| 



.4 



12 



'12 



V 1 - 9l2 



bt. 



A 



12 



12 



V 1 - 912 



with q u = qi + q2, leads to bosonic commutation rules for the operators bu, bf 2 - In 
analogy to the TRPA, we introduce phonon creation operators as 



Q+ = E^W 2 -0^i2 



(9) 



12 



and define the thermal vacuum state as the vacuum for the new phonon operators. 
Then 

|* (r)> = -L«p5|o(T)> 



i 



N 



S — - E Ci234&l2^34 
/ 1234 



where the matrix C1234 is defined through: 

^ ^12^1234 



y 31 



12 



For practical applications it is more appropriate to write the system of equations of 
the renormalized TRPA in terms of new variables, namely 



^24 



1 



^2 2/4 + 2/2^4 



(^24 - 042 ) 



^24 



X21J4 ~ V2%4 



12/ 



The result is 



E 



*12 + E ^1324^26^46 (^595 + ylO~ ~ Qs)) 



345 



26 



16 



(10) 



E ^12(^2 + ^(1- 92)) = jv 



11^ 



12 



^24*24 + ~ E E (V 567 8D* 51 D* 2 D 73 D 84 ) yjl - q 13y /l - q 24 ( ni - n 3 )^3i = 

Z 13 5678 



24 



9i = E [(^ i2^i) 2 + (%iV2Xi) 



v>0 



(12) 
(13) 



This system of equations is nonlinear. The important feature is that the mean-field 
variables not only dependent on T but also on the collective variables $. This 
coupling disappears in TRPA. If qi—0, the system splits into two independent parts: 
eqs. (10)-(11) and eq. (12). The latter then corresponds to the TRPA equation. 

Only the harmonic part of the thermal Hamiltonian is taken into account within 
the TRPA and renormalized TRPA. To incorporate effects of 7^31 and 7ii3 one has 



to extend the class of variational functions by including thermal four-quasiparticle 
configurations ||1 1|| : 



12 1>2,3>4 

Direct application of the variational principle to (14) yields the thermal second RPA 
matrix equation: 

\B* A * ) ( -Y ) = Uv ( Y ) • 
The stability matrix is composed of 2 x 2-supermatrices: 

A _ I ^12;34 ^12;3456 ] . B = ( ^ 12;34 ^12;3456 ] 

\ ^4l234;56 ^1234;5678 / \ -Bl234;56 -Bl234;5678 / 

and the amplitudes X u and Y v are two-component supervectors: 
X v = | ^ 2 | ■ Y v —[ 12 | 

\ ^1234 / V 1234 / 

Explicitly, the supermatrices A and B are expressed as ^(^-expectation values of 
double commutators of the thermal Hamiltonian Ji with different kinds of two and 
four thermal quasiparticle operators. 

If one approximates l^oCO) by |0(T)) (as in the TRPA), the matrix elements 
A 12; 34 and -E?i2;34 are the usual TRPA matrix elements while A 12; 3456, ^4i234 ; 56, an d 
^4i234;5678 describe the coupling of two-quasiparticle states with more complex ones: 

-4l2;34 = (£3 - ^4)^23^14 + ^2413 (^42/22/3 - x 2 x 3 y 1 y 4 ) 
-4l2;3456 = a(65)[/ 34 26 (^42/22/6 ~ ^62/22/32/4)5l5 + 

+a(3A)U 5613 (x 5 x 6 y 1 y 3 - x 3 y 1 y 5 y 6 )5 2 4 

^1234;5678 = ^45^36^27^18 ( £ 3 + e i ~~ £ 1 ~ £2) + 

+527518^4356(^3^43:5^6 ~ 2/3 2/4?/5 2/6 ) ~ 
-545^36^2178 {xiX 2 X 7 X 8 - y t J/22/72/8) + 

+a(21)a(78)a(43)a(56)5i 8 5 36 f/472 5 (a;4X5y22/7 - x 2 x 7 y A y 5 ) 

5l2;34 = U 3 214(x2X 3 y 1 y 4 - XiX^ys) 

where a(ij) denotes the antisymmetrizer, i.e. a(ij)Fij = Fij — F^. The coefficients 
-Bi234;56 and 51234J5678 vanish. 

In its full form the TSRPA is quite complicated for practical calculations. It is 
possible to simplify the problem by using of the boson operators A + ,A (b + ,b) or 



the phonon operators Q + ,Q constructed from A + ,A (b + ,b) as building blocks in 



the trial wave function (14) and the thermal Hamiltonian (6), as was done in [12" 



extending the quasiparticle-phonon nuclear model |13| to T 7^ 
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